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Abstract. We discuss some phenomenological aspects of y-ray emitting jets. In particular, we 
present calculations of the y-sphere and ^-sphere for various target photon fields and employ 
them to demonstrate how y-ray observations at very high energies can be used to constraint the 
Doppler factor of the emitting plasma and the production of VHE neutrinos. We also consider some 
implications of the rapid TeV variability observed in M87 and the TeV blazars, and propose a model 
for the very rapid TeV flares observed with HESS and MAGIC in some blazars, that accommodates 
the relatively small Doppler factors inferred from radio observations. Finally, we briefly discuss the 
prospects for detecting VHE neutrinos from relativistic jets. 
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INTRODUCTION 

Ejection of collimated relativistic outflows appears to be a common phenomena in as- 
trophysics. The radio-to- y-ray continuum emission observed in blazars, microquasars, 
and y-ray bursts (GRBs) is believed to be produced in such outflows on various scales. 
The common view is that these flows are powered by a magnetized accretion disk and a 
spinning black hole, and collimated by magnetic fields and/or the medium surrounding 
the jet. However, there is as yet no universal agreement about the mechanisms responsi- 
ble for the formation and acceleration the jet and the dissipation of its bulk energy. Even 
the composition of the jet is unknown in most sources. 

Active states during which rapid, large amplitude variations of the high-energy emis- 
sion are observed appear to be quite common in most classes of compact relativistic sys- 
tems. This activity is presumably associated with violent ejection episodes, as directly 
inferred in a few cases. Some of the objects mentioned above may also provide sites for 
acceleration of the UHE cosmic rays detected by various experiments, if the latter are 
indeed produced in a bottom-up scenario, posing a great challenge to the theory of parti- 
cle acceleration. Emission of VHE neutrinos should accompany the production of those 
UHECRs, and optimistic models predict fluxes in excess of detection limit of upcoming 
cubic-km scale neutrino telescopes. As discussed below, observations of VHE y-rays 
can provide stringent constraints on the photopion opacity, which can be translated into 
upper limits on the neutrino flux. 

Of particular interest is the class of TeV sources. There are at present over a dozen 
TeV blazars (for updated list see, e.g., Ref ill]), all of which are exclusively associated 
with the class of high peak BL Lac objects, another (nonblazar) TeV AGN, M87 [2], and 
several X-ray binaries (microquasars or gamma-ray binaries). The observed bolometric 
luminosity of TeV blazars during quiescent states is typically of the order of a few times 



10 ergs s , with about 10 percents emitted as VHE y-rays. The luminosity in the VHE 
band may be larger by a factor of 10 to 100 during flaring states. The intrinsic spectra 
(corrected for absorption on the extragalactic background light) appear to be hard, with 
a peak photon energy in excess of 10 TeV in the most extreme cases. The constraints on 
the dynamics of the system are most stringent in this class of sources, and are discussed 
in some greater detail below. 

STRUCTURE AND DYNAMICS OF y-RAY JETS 

The y-sphere and the 7T-sphere 

The pair production opacity is typically large within the inner jet region in essentially 
all classes of compact high-energy sources. This implies that y-rays produced at small 
radii will not be able to escape the system before being converted to e ± pairs. Both 
the synchrotron photons produced inside the jet and ambient radiation intercepted by 
the jet contribute an opacity to pair production. Results of detailed calculations of the 
pair production opacity are exhibited in figure 1, where the y-spheric radius, defined 
as the radius fy(£y) beyond which the pair production optical depth to infinity is unity, 
viz., T ry (r r ,£ r ) = 1, is plotted against y-ray energy e y , for two target radiation fields: 
synchrotron radiation (dashed lines) and external radiation (solid lines). The spectra 
of the target radiation fields employed in those calculations are given in Ref Q3[| . As 
seen the y-spheric radius increases, quite generally, with increasing y-ray energy, and 
for luminous sources can be much larger than the dissipation radius (indicated in the 
figure). 

In the powerful blazars, like 3C279, and in microquasars the intensity of both external 
and synchrotron radiation is quite large, corresponding to the uppermost curves in fig. 
1. This then implies that the y-spheres at energies corresponding to the GLAST band 
should encompass a rather large range of radii. In these sources the y-spheres can be 
mapped in principle by measuring temporal variations of the y-ray flux in different 
energy bands during a flare. If the y-ray emission is produced over many octaves of 
jet radius, where intense pair cascades at the observed energies are important 0], then 
it is expected that a flare will propagate from low to high energies, or that the variations 
at higher y-ray energy will be slower than at lower energies. With the limited sensitivity 
and energy band of the EGRET instrument it was practically impossible to resolve such 
effects. It is hoped that with the upcoming GLAST instrument this will be feasible. In 
contrast, in TeV blazars, which are much fainter, the target photon luminosity is small 
and the location of the y-spheres is not constrained, except, perhaps, at the highest 
energies observed (a few TeV). This difference should be reflected in the variability 
pattern of the VHE emission. A particular model for the rapid TeV flares is discussed 
below. 

The same target radiation field contributes also an opacity to photopion production. 
Because both the protons and the y rays interact locally with the same target radiation 
field the ratio of photomeson and pair production opacities depend solely on the ratio of 
cross sections, o p7 / a rr ~4x 10~ 3 , and the spectrum of the target radiation field. For a 
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FIGURE 1. Dimensionless y-spheric radius versus y-ray energy, computed in Ref J^t]. The different 
curves correspond to a different normalization of the target radiation field intensity. The dissipation radius 
rj is indicated. 
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Detailed calculations of opacity ratios that employed more realistic target photon spectra 
are presented in Ref [3]. It is found that at y-ray energies above a few TeV the opacity 
ratio is smaller than unity even at the maximum proton energy (determined from the 
confinement limit). For the class of TeV blazars this implies neutrino yields well below 
detection limit. With GLAST it should be possible to constrain other sources and to use 
such constraints to identify the best candidates for the upcoming km 3 detectors. 



Implications of variability 

The observed variability limits the linear size of the emission region (as measured in 
the Lab frame) to d ^ 10 14 (1 — /3 cos 0«) ? V arh/(l +z) cm, where /3 is the bulk speed 
of the emitting fluid, n is the viewing angle, ? var h is the observed variability time in 



hours, and z is the redshift of the source. The rapid variability observed in GRBs and 
blazars implies typically d/r g <T 2 , where T is the bulk Lorentz factor of the emitting 
fluid. The rapid VHE flare recorded recently in the TeV blazar PKS 2155-304 requires 
r ~ 20 in order that d ~ r g , regardless of any other considerations. 

The location of the emitting plasma is yet another issue. If the emission originates 
from radii r em ~ d, as often assumed, then compactness arguments yield a lower limit 
on the Doppler factor of the emission zone, as discussed further below. If, on the other 
hand, the emission is produced at radii r em >> d, as proposed recently for M87 and 
TeV blazars, then the fraction of jet energy that can be dissipated and converted to 
radiation in a conical jet of opening angle is at most r\ ^ (d/0r em ) 2 . As a consequence, 
either the opening angle of the jet must be very small, 0~J/r<<l,or the radiative 
efficiency L rac j/L ; must be very small, implying unreasonably large jet power in the most 
extreme cases. In situations where the jet is underpressured relative to the confining 
medium the jet is expected to converge to the axis. Reflection of collimation shocks 
at the nozzle may then give rise to appreciable dissipation in a very small region. The 
pattern speed of the emission region can differ significantly from the speed of the fluid, 
and stationary features, as occasionally observed in radio jets of blazars pD can be 
naturally produced. Such a model has been proposed to explain the rapid variability 
of the resolved X-ray emission and the unresolved TeV emission from the HST1 knot 
in M87 |@]. Alternative explanations have been offered for the rapidly varying TeV 
emission (e.g., Ref fl7J], and references therein). It should be noted though that in M87 
the X-ray and TeV luminosities, Lj e y ~ L x ^ 10 41 erg s _1 [SB], are much smaller 
than the TeV luminosity, Lx e v ~ 10 44 ~ 45 erg s _1 , observed typically in the TeV blazars. 
Estimates of the jet power in M87 yield Lj ^ 10 44 erg s _1 [8J], implying a very small 
conversion fraction, Lx e v /Lj ~ 10~ 3 . Even with such a small conversion efficiency an 
opening angle 6 < 10~ 2 rad is required if the TeV emission were to originate from 
the HST1 knot, unless reconfinement can give rise to sufficient convergence of the jet 
at the location of HST1, as proposed in Ref [6]. This idea is compelling since even 
modest radiative cooling of the shocked jet layer in a proton dominated jet will lead 
to such a convergence, at least in the non-relativistic case [9]. The effect of cooling on 
the collimation of relativistic jets needs to be explored. The stationary radio features 
observed in blazars seem to indicate that recollimation shocks may be an important 
dissipation channel in blazsrs, and this may apply also to other sources, e.g., GRBs 
[fioh . Whether the extreme TeV flares observed in VHE blazars can be accounted for by 
recollimation shocks at radii r em » d remains to be investigated. However, this would 
not resolve the 'Doppler factor crises' if the IR emission would turn out to vary on 
timescales comparable to the duration of the TeV flare. 



Constraints on Doppler factors 

Constraints on the Doppler factor of the fluid emitting VHE y-rays can be derived 
by measuring the low-energy flux (radio-to-IR) simultaneously with the variable y-ray 
emission. The requirement that the pair production opacity should not exceed unity, viz., 
T^°(r em , £y/ 8) < 1, constrains the density of target photons: n s (r em ) < (Oyyd)~ l . If the 



emission is assumed to originate from the innermost jet radii, in which case r em ~ d, 
then we must have ry(gy) < r em . The latter condition on r r (e r ) can be solved for the 
Doppler factor to yield yQ 

8 5 > 2 x 10 11 r; ar 3/ h 2 (re)-V(e r /lTeV) 1 /2 5jy , (2) 

where 5j y is the measured synchrotron flux density in Janskys and z is the redshift of the 
source. In deriving eq. © a luminosity distance di = 10 28 z has been adopted. In cases 
where r em » d the compactness of the TeV emission zone may be constrained by the 
variability of the IR flux observed simultaneously with the TeV flare. If the IR emission 
varies over time scales comparable to the duration of the VHE flare then eq. © still 
applies. Lower values of 8 are allowed if the variability time of the IR emission is much 
longer than the variability time of the VHE emission. 

Adopting TO ~ 1, we estimate 8 > 35 for the rapid flare observed in Mrk 421 and 
8 > 190 for the few minuets variability reported for PKS 2155-304, with r em < 10 17 cm 
for the minimum condition in both sources. 



A model for rapid flares in TeV blazars 

The large values of the Doppler factor implied by opacity constraints and rapid 
variability in TeV blazars are consistent with those obtained from fits of the SED to 
a homogeneous SSC model, but are in clear disagreement with the much lower values 
inferred from unification models llll [l2ll and superluminal motions on parsec scales 
tfULSHl. Various explanations, including a structure consi sting of interacting spine and 
sheath H$\ . opening angle effects iflot l and jet deceleration fl 1 Vt. 1 1 8h have been proposed 
in order to resolve this discrepancy. 

It has been proposed recently |l9|] that the rapid TeV flares observed in sources like 
Mrk 421, Mrk 501 and PKS 2155-304 are produced by radiative deceleration of fluid 
shells expelled during violent ejection episodes. These shells are envisaged to accelerate 
to a Lorentz factor To >> 1 at some radius r c \ ~ 10 2 — 10 3 r g , at which dissipation of their 
bulk energy occurs. The dissipation may be accomplished through formation of internal 
shocks in a hydrodynamic jet or dissipation of magnetic energy in a Poynting flux 



dominated jet [23,l2ll], and it is assumed that a fraction <^ e of the total proper jet energy 
density, uL is tapped for acceleration of electrons to a maximum energy y max w e c 2 . The 
dynamics of the front is then governed by the equation fl9ll 

!>;.r 2 /3) = -^y 4 ^, (3) 

where u s is the energy density of the target radiation field, as measured in the Lab 
frame, and % =< y 2 > /(< y > y ma x) depends on the energy distribution of nonthermal 
electrons. For a power law distribution, dn e /dy°^ y~ q with q < 2, we have 1 > % > 0. 1. 
Under the assumptions that u s (r) r~ 2 and that the proper density and average energy 
of the nonthermal electrons are independent of radius the solution of eq. © (in the limit 



/3 = 1) reads: 



^oy— , (4) 



where is the asymptotic Lorentz factor downstream. The stopping length can be 

By] 



expressed in terms of the optical depth for 77 absorption of a 7-ray of energy m e c 2 £y by 



a power law target photon field of the form I s (e s ) 

with g(e r ) = (CyC^) ^ 1 if a > 1 and g(e r ) = {e y e s ^ max ) a ^ 1 if a < 1, and g(e y ) < 
1 in both cases. We conclude that for a reasonably flat distribution of nonthermal 
electrons, q < 2, extension of the distribution to a maximum energy y max at which 
the pair production optical depth, T(To7 m ax), is a few is already sufficient to cause 
appreciable deceleration of the front. 

From the above it can be shown [19] that for the TeV blazars a background luminosity 
of L s ~ 10 41 — 10 42 erg s , roughly the luminosity of LLAGN, would lead to a substan- 
tial deceleration of the front and still be transparent enough to allow the TeV 7-rays 
produced by Compton scattering of the background photons to escape the system. The 
ambient radiation field is most likely associated with the nuclear continuum source. The 
bulk Lorentz factor of the jet during states of low activity may be appreciably smaller 
than that of fronts expelled during violent ejection episodes. 

WHAT IS THE PROSPECT FOR DETECTION OF VHE 
NEUTRINOS FROM RELATIVISTIC JETS? 

A new generation of experiments just started operating or will become operative soon, 
design to detect VHE neutrinos (IceCube, ANTARES, NESTOR, NEMO), and UHE 
cosmic rays (HiRes and the hybrid Auger detectors), are probing and will probe regions 
opaque to electromagnetic radiation, which are presently unaccessible, and will deter- 
mine the composition of jets. Besides providing an important probe of the innermost 
regions of compact astrophysical systems, these experiments can also be exploited to 
test new physics. 

As mentioned above, the detection of UHECRs motivated considerations of heavy 
jets that effectively accelerate protons to energies approaching the confinement limit. 
Effective neutrino production requires large photopion opacity which, as discussed 
above, can be constrained by VHE 7-ray observations. Present observations rules out 
TeV blazars as potential candidates for the upcoming km 3 neutrino telescopes (see 
Ref [[l] for further discussion). This leaves GRBs, microquasars, and perhaps powerful 
EGRET blazars as potential candidates. Powerful blazars at a redshift of z = 1 may 
produce up to one event per year in a km 3 detector [22] . In the case of microquasars 
up to a few events can be detected during a strong outburst if the viewing angle is 



sufficiently small [23,|24|,|25|]. The estimated neutrino flux from GRBs implies that only 



nearby sources can be individually detected by the upcoming experiments. However, the 



cumulative flux produced by the entire GRB population should be detectable assuming 
that cosmological GEBs are the sources of the observed UHECRs Il26ll . 
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